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Abstract. Lekkerkerker and Boland characterized the minimal forbid- 
den induced subgraphs for the class of interval graphs. We give a linear- 
time algorithm to find one in any graph that is not an interval graph. 
Tucker characterized the minimal forbidden submatrices of matrices that 
do not have the consecutive-ones property. We give a linear-time algo- 
rithm to find one in any matrix that does not have the consecutive-ones 
property. 



1 Introduction 

A graph is an interval graph if it is the intersection graph of a set of intervals on a 
line. Such a set of intervals is known as an interval model of the graph. They are 
an important subclass of perfect graphs [3], they have been written extensively 
about and they model constraints in various combinatorial optimization and 
decision problems |ll[|13j . They have a rich structure and history, and interesting 
relationships to other graph classes. For a survey, see [T]. 

If M is a 0-1 (binary) matrix, we let size{M) denote the number of rows, 
columns and I's. Such a matrix has the consecutive- ones property if there exists 
a reordering of its columns such that, in every row, the I's are consecutive. A 
clique matrix of a graph is a matrix that has a row for each vertex, a column for 
each clique, and a 1 in row z, column j if vertex i is contained in clique j . A graph 
is an interval graph if and only if its clique matrices have the consecutive-ones 
property, see, for example (2- 

In 1962, Lekkerkerker and Boland described the minimal induced forbidden 
subgraphs for the class of interval graphs [6 , known as the LB graphs (Fig- 
ure [5]). Ten years later. Tucker described the minimum forbidden submatrices 
for consecutive-ones matrices |16) . These are depicted in Figure [T] Not surpris- 
ingly, there is a relationship between the intersection graphs of rows of Tucker 
matrices and the LB graphs, depicted in Figure [H 

In this paper, we give a linear time bound for finding one of the LB subgraphs 
when a graph is not an interval graph. As part of our algorithm, we also give 
a linear-time {0{size{M)) bound for finding one of Tucker's submatrices in a 
matrix M that does not have the consecutive-ones property. This latter problem 
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Fig. 1. The Minimal Forbidden Submatrices for Consecutive-Ones Matrices. For Mi, 
fc > 3, and for Mu and Af/j, k > 4. Mjv and Mj have fixed size. 



was solved previously in 0{n * size{M)) time in where n is the number of 
rows of the matrix. 

A graph is chordal if it has no chordless cycle {Gm on four or more vertices). 
Figure [5] seems to imply that the rows of a given Tucker submatrix of the clique 
matrix can be extended to an interval model of an LB subgraph by including at 
most three additional rows of the matrix, giving the rows of vertices that induce 
an LB subgraph. Figure [5] gives a counterexample. Fortunately, it is true in the 
case of chordal graphs, but the example illustrates the need for a proof; it does 
not follow from seemingly obvious considerations, such as that no clique is a 
subset of any other. This gives an LB subgraph if G is chordal, and when it is 
not chordal, the algorithm of [15] already gives a Gm. 

An interval graph is proper if there exists an interval representation where no 
interval is a subset of another. It is a unit interval graph if there exists an interval 
representation where all intervals have the same length. These graph classes are 
the same, and Wegner showed that a graph is a proper interval graph if and only 
if it does not have a chordless cycle, the special case of G/y or Gy for n = 6 or 
the claw {Ki^^) as an induced subgraph [T7]. Hell and Huang give an algorithm 
that produces one of them in linear time [4] . The problem of finding a forbidden 
subgraph reduces easily to finding an LB subgraph. Each of the LB graphs 
is either one of Wegner's forbidden subgraphs or contains an obvious claw, and 
finding a claw in linear time, given an interval model, is elementary. By itself, this 
approach has no obvious advantages over Hell and Huang's elegant algorithm, 
but such reductions are useful when studying or programming a collection of 
related algorithms. 

A certifying algorithm is an algorithm that provides, with each output, a 
simple-to-check proof that it has answered correctly [5l |9]. An interval model 
gives a certificate that a graph is an interval graph, and an LB subgraph gives 
one if the graph is not an interval graph. However, a certifying algorithm was 
given previously in j5] . The ability to give a consecutive-ones ordering or a Tucker 
submatrix in linear time gives a linear-time certifying algorithm for consecutive- 
ones matrices, but one was given previously in [S]. The previous certificates are 
easier to check, which is a desirable property for certifying algorithms. However, 
they are neither minimal nor uniquely characterized. Aside from the theoretical 
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interest in LB subgraphs, it is easy to obtain a minimal certificate of the form 
given in [5] from an LB subgraph found by the algorithm we describe below. 
An open question is whether a minimal, unique, especially simple, or otherwise 
interesting special case of the certificate from [8] can be obtained by applying the 
algorithm of that paper to a Tucker submatrix obtained by the algorithm of the 
present paper. Tucker submatrices may be useful in heuristics for finding large 
submatrices that have the consecutive-ones property, small Tucker matrices, or 
identifying errors in biological data [Mj [2] . Our techniques provide new tools for 
such heuristics. 




A 



M,, 



n>3 
^ R,, 



k>3 




Ro R, R2 Rk-2 
k>2 



k>3 



Fig. 2. Gi through Gv are the minimal non-interval graphs discovered by Lekkerk- 
erker and Boland. Below them are the intersection graphs of the corresponding Tucker 
matrices. 



2 Preliminaries 

Given a graph G, let V denote the number of vertices and E denote the number 
of edges. If C X C F, let G[X] denote the subgraph induced by X. Standard 
sparse representations of 0-1 matrices take 0{size{M)) space to represent M. 
We treat the rows and columns as sets, where each row R is the set of columns 
where the row has a 1 and each column C is the set of rows where the column 
has a 1. Suppose TZ is the set of rows of a consecutive-ones ordered matrix and 
(Ci, C2, . . . , Cm) is the ordering of the columns. In linear time, we can find, for 
each row, the leftmost and rightmost column in the row. Let us call these the 
left endpoint and right endpoint of the row. 

That interval graphs are a subclass of the class of chordal graphs follows 
from inclusion of G777 among the LB subgraphs. Rose, Tarjan and Lueker give 
an 0{V + E) algorithm that finds whether a graph is a chordal graph, and, if 
so, produces its maximal cliques [12]. Otherwise, the algorithm of [T5^ produces 
a chordless cycle (G777) in linear time. 
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When a graph is chordal, the problem of deciding whether it is an interval 
graph reduces to the problem of deciding whether its clique matrix has the 
consecutive-ones property. Booth and Lueker further reduced this problem to 
that of finding a maximal prefix TV — {Ri, ■ ■ , Rr} of the rows of a binary 
matrix M that has the consecutive-ones property, in 0{size{M)) time. 

Assigning a left-to-right order to children of each internal node of a rooted 
tree results in a unique left-to-right order of the leaves. Booth and Lueker's algo- 
rithm produces a PQ tree, for TZ' . The PQ tree represents all possible consecutive- 
ones orderings of TZ' . There is one leaf {c} for each column c. The internal nodes 
of the PQ tree consists of P nodes and Q nodes. The consecutive-ones ordering 
of columns are given by the leaf orders obtainable by assigning an arbitrary left- 
to-right order to children of each P node, and for each Q node, assigning the 
given left-to-right order or its reverse. 

Though the PQ tree can be represented using 0(1) space per node, concep- 
tually, we will consider each node of the PQ tree to be set given by the disjoint 
union of its children; equivalently, it is the union of its leaf descendants. 

Definition 1. Let S be a collection of subsets of a set U . Two elements of U 
are in the same Venn class if they are elements of the same set of members of 
S. The unconstrained Venn class consists of those elements of U that are not in 
any member of S; all others are constrained. Two sets Ri,R2 overlap if their 
intersection is nonempty, but neither is a subset of the other. The overlap graph 
ofS is the undirected graph whose vertices are the members of S , and Ri, R2 € S 
are adjacent if and only if Ri and R2 overlap. 

Lemma 1. ^10] A set of columns is a Q node of a consecutive- ones matrix M 
if and only if it is the union of rows of a connected component of the overlap 
graph of rows of M . The Venn classes of rows in this component are its children. 

3 Breadth-first search on the overlap graph of a collection 
of sets, given a consecutive-ones ordering 

In linear time, we may label each row of a consecutive-ones ordered matrix with 
its left and right endpoints. We may then label each column Ci of a consecutive- 
ones ordered matrix with the set of rows that have their left endpoints in c^. In 
linear time, we can then radix sort the list of sets that have their left endpoint 
at Ci in descending order of index of right endpoint, yielding a list TZi. This is 
accomplished with a single radix sort that has the index of the left endpoint as 
its primary sort key and index of the right endpoint as the secondary sort key. 
By symmetry, we can construct a list Ci of rows that have their right endpoint 
in each column Ci, sorted in ascending order of index of left endpoint. 

This allows us to perform a breadth-first search on the overlap graph of the 
rows in time linear in the size of the matrix, as follows. The lists Ci and TZi are 
represented with doubly-linked lists. We maintain the invariant that elements 
that have been placed in the BPS queue have been removed from these lists. 
When a consecutive-ones ordered set R comes to the front of the queue, we 
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traverse its list (cj,Cj+i, . . . ,Ck) of columns. For each Ch in the list, we remove 
elements from TZ^ and place them in the queue, until we reach an element in 
TZi whose right endpoint is no farther to the right than Cfc. All of the removed 
elements overlap R. Since TZh is sorted in descending order of right endpoint, 
all these elements are a prefix of TZh , and any remaining elements in the list do 
not overlap R. When we remove an element from TZh, we remove it from any 
list Ch' that it is a member of, to maintain the invariant. This takes 0(1) time 
for each element moved to the BFS queue, plus 0(1) time for each column of R. 
The lists Lh are handled symmetrically. Summing over all rows R, the time is 
0{size{M)). 

4 Finding Tucker Submatrices 

4.1 Tucker matrices with at most four rows 

Lemma 2. // a set TZ' of rows has the consecutive-ones property and Z is a row 
such that TZ = TZ' U {Z} does not, then Z is one of the rows of every Tucker 
submatrix in TZ. 

Algorithm 1 (See Lemma\^) 

initialRows(M' , k ) 
M M' 
I ■= 1; 

While i < fc + 1 and M has at least i — 1 rows 

Using Booth and Lueker, find the minimal prefix R2, . . . , Rr, Z) of rows 
of M that does not have the consecutive-ones property. 

Let M be the matrix whose row sequence is {Z, Ri, R2., . . . , Rr). 
return M 

Lemma 3. Suppose Algorithm[l\is run with parameter k and a matrix M' that 
does not have the consecutive-ones property. If the returned matrix M has at 
most k rows, then these are the rows of every Tucker matrix of M . Otherwise, 
M fails to have the consecutive-ones property and every Tucker submatrix in M 
has at least fc + 1 rows. 

Proof. By induction on i, M does not have the consecutive-ones property at 
the end of iteration i. Also, by induction on i, using Lemma [21 at the end of 
iteration j, for every Tucker submatrix Mt of M, the rows of Mt include the 
first i rows of M. If Af^ has only i rows, then the first i rows of M do not have 
the consecutive-ones property, so at the end of iteration i + \, M will have i 
rows. 

We run Algorithm [T] for fc = 4. If it returns a matrix with j rows, where j < 4, 
it is easy to get a linear time bound to get the columns. (One way is to generate 
all j! < 24 orderings of rows and for each, to check for the columns of each 
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Tucker matrix of size j.) Otherwise, Algorithm [T] returns a matrix M of more 
than 4 rows. By Lemma [U M fails to have the consecutive-ones property and 
every Tucker submatrix of M has at least five rows. This excludes any instances 
of M/y, My or the anomalous case of M/ on three rows that does not correspond 
to a chordless cycle. 

4.2 Matrices in which all Tucker submatrices have more than four 
rows 

Lemma 4. The overlap graphs of Mi, Mu, and Mjjj are simple cycles. 

Definition 2. Suppose TV is a set of rows with the consecutive- one property, Q 
is a Q node of its PQ tree, {Xi,X2, ■ ■ ■ , Xk) is the ordering of Q 's children and 
Z is a row not in TV . Let Xh,Xi,Xj be three children of Q such that h < i < j. 
They are a 1-0-1 configuration for Z if Xh and Xj each contain a 1 of row Z 
and Xi contains a of row Z. They are a 0-1-0 configuration for Z if X^ and 
Xj each contain a of row Z and Xi contains a 1 of row Z . 

Lemma 5. IfTZ' is a set of rows that has the consecutive-ones property and Z 
is a row not in TV , then TV U {Z} does not have the consecutive-ones property if 
the PQ tree of TV has a Q node Q such that either: 

1. Q has a 1-0-1 configuration for Z ; 

2. Q has a 0-1-0 configuration for Z and Z is not a subset of Q. 

This test is implicit in Booth and Lueker's algorithm, where it is a sufficient 
condition, but not a necessary one. The following is a consequence of Lemma [TJ 

Lemma 6. JS^ The conditions of Lemma\5\ are necessary and sufficient if the 
overlap graph of of TV is connected. 

Lemma 7. // a matrix fails to have the consecutive- ones property and has no 
Tucker submatrix with fewer than five rows, then when Algorithm]^ is run on it 
with k — A, at the end of one of the five iterations of its loop, the PQ tree of 
TZ' ~ {-Ri, -R2, ■ • ■ , ^r} will have a Q node Q with the following properties: 

- Q has a 1-0-1 configuration (Xh, Xi, Xj) for Z; 

- There exist A, B € TZ' that are members of the component of the overlap 
graph on TV whose union is Q, and such that A contains Xh and is disjoint 
from Xi and Xj , and B contains Xj and is disjoint from Xh and Xi . 

Proof. If T is the rows of A/ that contain a Tucker submatrix Af^, at the end 
of an iteration of the loop of Algorithm [U Z £ T by by Lemma [21 By Lemma lU 
the overlap graph of T' = 7"\ {Z} is connected, so T' is a subset of a component 
of the overlap graph of TZ' , which gives rise to a Q node Q of the PQ tree of TV , 
by Lemma [TJ Since the children of Q are the Venn classes of the component, no 
two Venn classes of T' , hence no two columns of Af^, can lie in the same child 
of the Q node. 
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For each choice of a last row of a Tucker matrix on at least five rows, Fig- 
ure |3] gives the possible orderings imposed on the last row by a consecutive-ones 
ordering of T', which is unique up to reversal, by Lemmas |4] and [TJ In each case, 
if row i ^ {0, 1, fc — 2, fc — 1} is chosen to go last, rows i — 1 and « -I- 1 satisfy 
the requirements of A and B. Mr has at least five rows and no row of Mt is 
contained in Z more than once in the five iterations, so in at least one of the 
iterations a row i ^ {0, 1, fc — 2, fc — 1} will go last. 
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Fig. 3. Consecutive-ones orderings of all but the last row of Mj, Mu and Mm for 
different choices of last the last row. For all but at most four choices of the last row, 
there exists a 1-0-1 configuration and rows A and B satisfying Lemma [7] 



The correctness and linear time bound for the following are the key results 
of this section: 

Algorithm 2 Find the rows of a Tucker submatrix when every Tucker subma- 
trix has at least five rows 

findRows(M ) 

Run initialRows(M,4,) (Algorithm]!^ to find TZ' , Z, A, B 

satisfying the requirements of Lemma 
Let P be a shortest path in the overlap graph of TZ' from A to B 
Let Pi be a minimal prefix of P such that the union of {Z} and the 

set Vi of rows of Pi does not have the consecutive- ones property. 
Let P2 be a minimal suffix of Pi such that the union of {Z} and the 

set 1^2 of rows of P2 does not have the consecutive-ones property. 
Return V2 U {Z} . 

Lemma 8. If M does not have the consecutive- ones property and every Tucker 
matrix of M has at least five rows, then Algorithmic returns the set of rows of 
a Tucker matrix of M . 
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Proof. Since A and B lie in the same component of the overlap graph of 7?.', P 
exists. Since TZ' has the consecutive-ones property, so does V. Because V has 
a connected overlap graph, P, IJ'P is a single Q node of the PQ tree of V, by 
Lemma [TJ Because of A and B, X^, Xi, and Xj are contained in distinct Venn 
classes of V, and the ones containing Xh and Xj are constrained. Since [JV 
is consecutive, it must have a row that contains the Xi, hence the Venn class 
of V containing Xi is also constrained. Therefore, P U {Z} does not have the 
consecutive-ones property by Lemma [5l and Pi and P2 exist. By Lemma [21 all 
Tucker matrices in TZ' U {Z} contain Z, so this applies also to P2 U {Z}. 

Suppose there is a proper subset TZ" of the rows on P2 such that TZ" U {Z} 
contains a Tucker matrix. The overlap graph of TZ" is connected, by Lemma |4l 
Since P2 is a shortest path, it is a chordless path, so TZ" is a subpath of P2 
by Lemma ID Let TZ[ be the rows on Pi , excluding the last row on Pi . Let TZ'2 
be the rows of P2, excluding the first row on Pj. By the minimality of Pi and 
P2, TZ'i U {Z} and TZ'2 U {^} have the consecutive-ones property. Since TZ" is a 
subpath of P2, TZ" C TZ'i or TZ" C TZ'2, ^" U {^} has the consecutive-ones 
property, contradicting our assumption that it does not. Therefore, P2 U {Z} is 
the set of rows of a Tucker matrix. 
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ones property 



Figure S] gives an example on which we illustrate some implementation de- 
tails. Z is given by the I's and O's above the column numbers in the figure on 
the left, and TZ' is depicted by the intervals. The rows labeled A and B satisfy 
the requirements of A and B for Lemma [3 and P = {A, E, F, G, H, J, L, B) is 
a shortest path from A to P in the overlap graph of TZ' , found using the BFS 
algorithm of Section [31 

Using Booth and Lueker's terminology, we maintain labels on each class indi- 
cating whether it is full (contains only I's of Z), empty (contains only O's of Z), or 
partial (contains both I's and O's of Z). The minimal prefix Pi of P whose rows, 
together with Z, do not have the consecutive-ones property, is (A, P, P, G, i?, J). 
This is detected as follows. It is easy to verify that its sequence of constrained 
Venn classes is ({0, 1},{2},{3}, {4,5}, {6 - 9}, {10, 11}, {12}, {13,14}, {15}, 
{16}), and their full/partial/empty labels {F, F, F, F, P, P, E, E, E, E), respec- 
tively. Selecting a 1 from a full class, a from the first partial class and a 1 from 
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the second partial class gives a 1-0-1 configuration satisfying Lemma O It is the 
minimal such prefix. A smaller prefix, {A, E, F) has a 0-1-0 configuration, but it 
does not satisfy Lemma [5] because Z C AU E U F. 

The minimal suffix P2 of Pi that satisfies Lemma [5] is {F, G, H, J) , which is 
found in the same way by working on the reverse of Pi. Its sequence of con- 
strained Venn classes are ({2 - 9}, {10, 11}, {12 - 14}, {15}, {16}), labeled 
{P,P,E,E,E), respectively. Selecting a from the first partial class, a 1 from 
the next, and a from an empty class gives a 0-1-0 configuration. It satisfies con- 
dition [2] of the lemma, because the unconstrained class, {0,1} is partial, hence 
Z <^FUGUHUJ. 

Therefore, {F, G, H, J, Z} is the set of rows of a Tucker submatrix. A minimal 
set of columns that illustrates that it satisfies the lemma is {0, 6, 10, 12, 15, 16}. 
On the righthand side of Figure |4] is the resulting Tucker matrix, which matches 
the final configuration in the sequence for Mm in Figure |3l 

This example shows that the key to finding Pi and P2 is maintaining the 
sequence of constrained Venn classes and their full/partial/empty labels as rows 
are added in the order in which they occur on P or on the reverse of Pi. Since 
they are added in an order such that every prefix of the order has a connected 
overlap graph, the sequence is uniquely constrained after each row is added, by 
Lemma [1] When a row Ri is added, if it overlaps a constrained Venn class X, 
X must be replaced in the sequence with two Venn classes, {X \ Ri, X D Ri) or 
with [X n Ri, X \ Ri), whichever is required to maintain consecutiveness of Ri. 
If Ri intersects the unconstrained class, S, then Ri D S must be added at one 
extreme end of the sequence, whichever maintains consecutiveness of Ri. Details 
are given in [8]. 

The difference between this algorithm and that of ^ (and Booth and Lueker) 
is that, instead of testing at each iteration whether the next row Ri can be added 
to those considered so far without undermining the consecutive-ones property, 
it must repeatedly perform this test on the fixed row Z after each row Ri is 
added. We already know that Ri can be added, since TZ' has the consecutive- 
ones property. Like Booth and Lueker, the previous algorithm of [5] applies the 
full/partial/empty labels for Ri to facilitate the test, in 0{\Ri\) time, and then 
removes them before considering the next row Ri+i. Though we must perform 
the test on the fixed row Z at each iteration, instead of on Ri, we must do it 
0{\Ri\) time, not 0(|Z|) time, in order to retain the linear time bound. To do 
this, we leave the full/partial/empty labelings for Z from one iteration to the 
next, so that we only have to update them, using Ri, rather than re-creating 
them each time a new row is considered. 

To facilitate this, we keep updated labels c{X) and n{X) on each Venn class 
X, where c{X) denotes the cardinality of X and n{X) is the number of elements 
of Z in X. Labels only need to be updated when a Venn class is split. It is split 
into XnRi and X\Ri. We may find c(Xni?i) and n{XriRi) by counting them 
directly, since there are 0{\Ri\) of these elements. The classes are implemented 
with doubly-linked lists, and these sets are removed from the list for X, leaving 
it to represent X \ Ri. Subtracting c{X n Ri) and n{X n Ri) from the old labels 
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c{X) and n{X) gives the updated labels for c{X \ Ri) and n{X \ Ri) in 0(1) 
time. Each of the new classes is full if its c() and n() labels are equal, empty if 
its n() label is 0, and partial otherwise. 

To evaluate whether one of the conditions of Lemma [S] holds, it is easy 
to see that it suffices to keep track of transition pairs, which are consecutive 
pairs such that one contains a and one contains a 1. This happens when their 
full/partial/empty labels are unequal, or else both partial. When a new tran- 
sition pair forms, we have touched at least one member of the pair within our 
0(|i?i|) operations, so keeping track of these does not affect this time bound. 

Since finding P takes linear time by the BFS of Section [3l it remains only 
to bound the time required to find the first step, finding the elements A and B 
of Lemma [71 This is much more straightforward, since we apply it once for each 
iteration of the loop of Algorithm [TJ hence we can afford to take 0{size{M)) 
time for the test. We can apply the entire set of full/partial/empty labels for Z 
to the PQ tree within this bound, by working from the leaves to the root. 

For each Q node Q, the members of the overlap component whose union is Q 
are unions of more than one and fewer than all of its children. How to find them 
in linear time for all Q nodes has been described previously, for example in [7]. 
We find the rows of the overlap component that contain a Venn child of Q that 
is labeled full or partial (a "1"). Out of all such rows, let A' be the one with a 
leftmost right endpoint, and let B' be the one with the rightmost left endpoint. 
By a simple greedy swapping argument, the overlap component giving rise to 
Q contains an A and a B satisfying Lemma [7] if and only if A' and B' satisfy 
it, which happens if and only if there is a child between the right endpoint of 
A' and the left endpoint of B' that is labeled partial or empty (a "0"). This 
can also clearly be implemented so that the time bound over all Q nodes takes 
0{size{M)) time. 

Once the set 7^2 U {Z} of rows of a Tucker matrix have been found, it remains 
to find the columns. This is a set of columns, the removal of any one of which 
would undermine the conditions of Lemma [Sj which are satisfied initially by 

U {Z}. Deletion of a column undermines the lemma if and only if does one of 
the following: 

1. It disconnects the path in overlap graph on TZ] 

2. It undermines the only remaining 1-0-1 configuration, or 0-1-0 configuration 
with a 1 in the unconstrained class. 

The second test is elementary and omitted because of space constraints. For 
the first test, recall that P2 = (i?i, -R2, ■ ■ ■ ,Rk) is a chordless path in the overlap 
graph. Let A ^ {Ri \ R2] U {X\X ^ R^ n R,+i for i £ {1, 2, . . . , /c - 1}} 
U{y |y — Ri+i \ i?i for i e {1, 2, . . . , fc — l}}. Each element of A is consecutive- 
ones ordered, the sum of cardinalities of sets in A is 0{size{M)). The overlap 
graph remains connected if and only if every member of A contains at least one 
retained column. We give each column a list of members of A it is contained in 
and keep a counter on each element of A indicating the number of remaining 
columns it contains. When removing a column C, the counters can be updated 
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by decrementing the counters of members of A in its list. A column cannot be 
removed if removing it would decrement a counter to 0. 

5 Finding a Lekkerkerker-Boland Subgraph 

Tucker observed that the smallest graphs whose clique matrices contain a Tucker 
matrix must be exactly the LB graphs \TW . However, making use of this to find 
an LB subgraph is not as straightforward as it appears. It is easy to believe from 
Figure [2] that that the rows of every Tucker submatrix in a clique matrix can 
be extended to a clique matrix of an LB subgraph by including at most three 
additional rows. The rows of the result would therefore identify the vertices that 
induce an LB subgraph. 

That this reasoning is flawed is illustrated by Figure [S] Neither does it follow 
from the fact that no column of the clique matrix is a subset of any other; we 
discovered the example in the figure when trying to prove it using only this 
fact. Therefore, though the conclusion is true in the case of chordal graphs, this 
requires proof. 



11 11 
1 1 1 

_L \] 1 

1 1 1 
1 1 1 
1 

1 



Fig. 5. A Tucker matrix. Mm, in the clique matrix of a graph G. The only LB subgraph 
is the chordless cycle (b, c, e, d), which does not contain row a of the Mm. 

If G is not chordal, we may return a Gm by the algorithm of (TS] . Henceforth, 
we may assume that the graph is chordal. A clique tree of a chordal graph is 
a tree that has one node for each maximal clique, and with the property that 
for each vertex v of G, the cliques that contain v induce a connected subtree. 
Every chordal graph has a clique tree, see for example [3] . A vertex is simplicial 
if its neighbors induce a complete subgraph. The following are immediate from 
results that appear there: 

Lemma 9. Let T be a clique tree for a chordal graph G, and let K he a leaf. 
Then K contains a simplicial vertex of G. Let S be the simplicial vertices of K, 
and let T' be the result of deleting leaf K from T. Deleting S from G yields an 
induced subgraph that has T' as a clique tree. 

Definition 3. By shrinking a clique tree T , let us denote the operation of delet- 
ing the set S of simplicial vertices in a leaf K of T , yielding a graph with the 
smaller clique tree described by the lemma. 
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Lemma 10. If G is a chordal graph, T is a clique tree ofG, and G[X] is con- 
nected, then the cliques of T that contain members of X induce a connected 
subtree of T . 

Lemma 11. Let T be a clique tree of a chordal graph G, JC a collection of cliques 
of G, and a clique C ^ K-. If every member of K. — C contains a vertex that C 
does not have and G'[IJ K.\C] is connected, then C does not lie on the path in 
T between any pair of cliques of JC — G . 

Definition 4. // a Tucker matrix occurs as a submatrix induced by column set 
C and row set TZ of a clique matrix, let a private row for G G C be a row R ^ TZ 
that is contained in no column of C other than G. Let the special columns of 
Mil — My be those that are subsets of other columns of these matrices, and let 
the special columns of an instance of Mj on three rows be all three columns. 

The special columns are the ones that Tucker identified as "asteroidal column 
triples" in the bipartite incidence graph of rows and and columns in [16]. How- 
ever, the existence of simplicial vertices or "asteroidal vertex triples" in a a sense 
defined in [6] was not examined. The paper dealt with arbitrary consecutive-ones 
matrices where such rows need not occur. 

Lemma 12. Let G be a chordal graph and let Mt be a submatrix of a clique 
matrix that is an instance of Mj on three vertices or an instance of A/// — My ■ 
Then the clique matrix has a private row for each special column of Mt- 

Proof. Let IC be the set of cliques of G that contain the columns of this instance 
of Mt, let C G /C contain a special column of the instance, and let T be a clique 
tree for G. In each case, G does not lie on the unique path in T between any 
pair of members of /C — C, by Lemma 1111 Therefore, iteratively shrinking T 
(Definition |3]) until it cannot be further shrunk without deleting a clique of K, 
results in an induced subgraph G' with a clique tree T' where C is a leaf. By 
Lemma [9l G has a simplicial vertex, which must be a private row for G in /C. 

By examination, it is easy to verify in the case of Mi on three rows, and 
Mil — Mjv tbat adding the private rows for the special columns to the rows 
and columns of the submatrix gives an interval model for the corresponding LB 
subgraph. The resulting set of rows therefore induces an LB subgraph. 
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